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We present a numerical study of a generalized two-dimensional Swift-Hohenberg 
model of spiral pattern formation in Rayleigh-Benard convection in a non-Boussinesq 
fluid. We demonstrate for the first time that a model for convective motion is able 
to predict in considerable dynamical detail the spontaneous formation of a rotating 
spiral state from an ordered hexagon state. Our results are in good agreement with 
recent experimental studies of C0 2 gas. The mean flow and non-Boussinesq effects 
are shown to be crucial in forming rotating spirals. 
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One of the most striking examples of spatio-temporal self-organized phenomena in 
non-equilibrium systems is the rotating spiral seen in chemical and biological systems 
HI] @. The Belousov-Zhabotinsky (BZ) reaction || has received considerable atten- 
tion as an example of chemical wave propagation. The formation of spiral patterns 
in the BZ system resulting from the coupling of a reaction process with a transport 
process such as diffusion has been extensively studied, both theoretically and experi- 
mentally, during the past decade. Remarkably, similar rotating spirals were observed 
recently in Rayleigh-Benard (RB) convection in non-Oberbeck-Boussinesq fluids [|J. 
Much of the earlier experimental work has been restricted to Oberbeck-Boussinesq 
type fluids, in which one observes various configurations of roll patterns. However, in 
a non-Oberbeck-Boussinesq system with, for example, temperature dependent trans- 
port coefficients, both roll and hexagonal patterns can exist. Very recently Boden- 
schatz et al 0] have performed experiments on convection in CO2 gas and studied 
the existence of and transitions between convective patterns exhibiting different sym- 
metries. They have observed the competition between a uniform conducting state, a 
convective state with hexagonal symmetry, and a convective state consisting of rolls. 
Their most surprising discovery is that the hexagon-roll transition has a tendency to 
form rotating spirals. In this paper we present the first numerical evidence for the 
spontaneous formation of a rotating spiral pattern during the hexagon-roll transition 
in a large aspect ratio cylindrical cell near onset. We also report the numerical re- 
sults of spiral formation during the conduction-roll transition. We will qualitatively 
compare our numerical results with the experiments of Bodenschatz et al. 

We now describle our two-dimensional (2D) study of spiral patterns formed in 
a thin layer of a non-Oberbeck-Boussinesq fluid. We model such a fluid by a two 
dimensional generalized Swift-Hohenberg equation || ||, given by equations (1) and 
(2) below, which we solve by numerical integration. The Swift-Hohenberg equation 
and various generalizations of it have proven to be quite successful in explaining many 
features of convective flow H, particularly near onset. As we show in this paper, the 
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same holds true for the non-Oberbeck-Boussinesq fluid. Our model is defined by 



dip(r, t) 
dt 



+ g m U ■ 



V 2 



ij) - g 2 i? + /(f), 



v 2 e = x v(v 2 v>) 



(2) 



where, 



[/ = {d v i)e x - {d x £)e y 



(3) 



with boundary conditions, 

M B = n- V^|b = 0, (4) 

where n is the unit normal to the boundary of the domain of integration, B. This 
equation with g 2 — g m — reduces to the Swift-Hohenberg (SH) equation and has 
been extensively used to model convection in thin cells and near onset |7], || |10| . 
The scalar order parameter ip{r, t) is related to the fluid temperature in the mid-plane 
of the convective cell. £{r,t) is the vertical vorticity potential. This mean flow field 
coupling II |ll| with the SH equation has been shown to play a key role in the onset 



of turbulence in Oberbeck-Boussinesq systems |12[ . The quantity e is the reduced 
Rayleigh number, 

e= — -l, (5) 

where R is the Rayleigh number and R™ is the critical Rayleigh number for an 
infinite system. A phenomeno logical forcing field / has been included in Eq.(l) to 
simulate the lateral sidewall forcing produced by horizontal temperature gradients 
present in the experiment. As in earlier studies |l3|, 14], we have varied the strength 



and spatial extent of / in order to best fit the experimental observations. We have 
derived a three mode amplitude equation from the generalized SH equation in order 
to both estimate the threshold values of e that separate regions in which roll and 
hexagonal configurations are stable, as well as the values of the parameters that enter 
the generalized SH equation in terms of experimentally measurable quantities. From 



the experiment [ffl, we find that #2 ~ 0.35, which is the value that we have used 
in Eq.(l). The value of e used in the numerical simulation is related to the real 
experimental value e exp in ref. 0] by e exp = 0.3594e. The nonlinear coupling constant 
g m has been chosen to be 50, which is consistent with earlier studies [JQJ and 
c 2 = 10 throughout our calculations. We note that Bestehorn et al ]15| have reported 
a study of this model limited to the special case in which the initial configuration is 
a spiral. This, however, avoids the fundamental question addressed here of how the 
spiral spontaneously forms in the hexagon-roll transition. In the following we report 
the results of our calculations. 



i) Formation of hexagonal pattern with sidewall forcing field. We have considered 
a circular cell of radius R = 32ir, which corresponds to an aspect ratio of T — 2R/n = 
64. A square grid with N 2 nodes has been used with spacing Ax = Ay = 647r/N, and 
iV = 256. We approximate the boundary conditions on ip by taking ip(f,t) = for 
||r|| > R, where r is the location of a node with respect to the center of the domain 
of integration. The initial condition ip(f,t = 0) is a random variable, gaussianly 
distributed with zero mean and a variance of 10 _1 . In this case e = 0.1, and / = 0.0 
everywhere except on the nodes adjacent to the boundary with / = 0.1. Figure 1 
presents a typical configuration for a hexagonal pattern which forms in the presence 
of a strong static sidewall forcing. 

ii) Early stage of transition between hexagons and rolls. We use Fig. 1 as the 
initial configuration with exactly the same forcing field / as before, but now we 
increase e very slowly up to 0.3. We take e = 0.1 + 1.67 x 10~ 4 t for < t < 1200 and 
e = 0.3 for t > 1200. Figure 2 shows two configurations during the early transient 
regime during the hexagon to roll transition. How the sidewall forcing and defects 
mediate the transition can be clearly seen in Fig. 2. The rolls are formed near 
the sidewall with a favorable orientation relative to the symmetry of the sidewall. 
In the meantime the defects glide toward each other and invade nearby regions of 
hexagonal order to create a region of rolls that spreads across the cell as the transition 
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proceeds. The spiral formation is already noticeable in Fig. 2(b). This resembles the 
experimental observation that there is a tendency to form spirals during the hexagon 
to roll transition. 



iii) Formation of rotating spirals. Fig. 3 show the spatial and temporal forma- 
tion of the spiral pattern at later times than in Fig. 2. In Fig. 3(a), we see that 
rolls bend or curve rapidly forming a roughly uniform patch of rolls with a locally 
disordered texture near the left corner. Further evolution consists primarily of dislo- 
cations gliding toward each other and eventually annihilating themselves, ending in 
a three- armed spiral. The final state of the rotating spiral (Figs. 3(b) and 3(c)) is 
remarkably similar to one observed in the experiment and occurs at t ~ 49000 ~ 12 
horizontal-diffusion times [PJ . The corresponding experimental times are in the range 
of 10 to 20 horizontal-diffusion times. Our numerical investigation indicates that the 
non-Boussinesq effect plays a crucial role in forming spontaneous spirals. In the ab- 
sence of $2 (with or without the mean flow field), if we start with a random initial 
condition, there is no occurence of a spiral pattern. This strongly suggests that the 
formation of spiral patterns is an intrinsic property of non-Boussinesq systems. We 
have studied the transition from a conduction state to a rotating spiral state, as well 
as the transition from a rotating spiral state to a hexagon state. We have also studied 
the effect of decreasing e in the hexagon to rotating spiral transition. We find that a 
stable n-armed spiral tends toward one with fewer arms when e is decreased, in agree- 
ment with experimental observation. The details of all the above will be discussed 
elsewhere. 



In summary, we have investigated the question of pattern formation in a model 
of convection in a non-Boussinesq fluid that allows patterns of various symmetries. 
We start with a random initial condition and show that this leads to the ordered 
hexagonal state observed in the experiment. We then show that upon increasing e we 
see a dynamical evolution to a new roll state which contains a rotating spiral pattern. 
These results are in very good agreement with the experimental studies in CO2 gas. 
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Figure captions 

Figure 1. Hexagonal pattern starting from random initial condition obtained in 
a cylindrical cell with aspect ratio T = 64. The values of the parameters used are 
g 2 = 0.35, g m = 50 and e = 0.1. A non-zero forcing field localized at the boundary 
with / = 0.1 has been used. 

Figure 2. We observe an early stage of hexagon-roll transition obtained by chang- 
ing e from e = 0.1 to e = 0.3, in a cylindrical cell with an aspect ratio T = 64. The 
initial condition is the uniform hexagonal pattern shown in Fig. 2. Two different 
times, t = 720 (b) and t = 960 (b) are shown. The rolls appearing near the defects 
and sidewall boundaries spreads through the cell as the transition proceeds. 

Figure 3. Pattern formation of a rotating three-armed spiral obtained in a cylin- 
drical cell with aspect ratio T = 64, with g 2 = 0.35, g rn = 50, e = 0.3 and f=0.1. 
Time series of the pattern evolution at (a) t= 7440, (b) t=48240 and (c) t=64080. 
The final rotating spiral pattern is shown in (b) and (c). 



9 



